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Abstract 

We study the deformed Harry Dym and Hunter-Zheng equations with two arbitrary deformation 
parameters. These reduce to various other known models in appropriate limits. We show that both 
these systems are bi-Hamiltonian with the same Hamiltonian structures. They are integrable and 
belong to the same hierarchy corresponding to positive and negative flows. We present the Lax 
pair description for both the systems and construct the conserved charges of negative order from 
the Lax operator. For the deformed Harry Dym equation, we construct the non-standard Lax 
representation for two special classes of values of the deformation parameters. In general, we argue 
that a non-standard description will involve a pseudo-differential operator of infinite order. 
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1 Introduction: 



The search for exactly solvable equations has acquired enormous importance since the nonlinear 
Kortweg-de Vries equation was shown to be integrable 0111 • Exactly solvable equations, linear or 
nonlinear, constitute a very special class of dynamical systems with many interesting properties. 
The significance of determining new integrable systems, as well as studying their properties and 
solutions, can not be overestimated. The structure of integrable systems (or partial differential 
equations) is highly restrictive and, in general does not allow for deformations. In particular, the 
presence of arbitrary constant parameters (which can not be transformed away by some symmetry) 
is quite rare. One notable exception is the two boson equation [3], which has an arbitrary constant 
parameter present, is known to be integrable. Different values of this parameter reduce the model 
to other known soluble models. In this paper, we propose a new system of equations, with two 
arbitrary constant parameters, that is exactly soluble and reduces to various other known physical 
models in different limits of these parameters. 

Let us recall that the Harry Dym (HD) equation ^ |5] 

= fc'/'). , (1) 

and the Hunter-Zheng (HZ) equation |S] 

{ut + uu^)^ = ^ul , (2) 

are known to be integrable and, in fact, belong to the same hierarchy corresponding to negative and 
positive order flows, respectively [7]. The HD equation has interest in the study of the Saffman- 
Taylor problem which describes the motion of a two-dimensional interface between a viscous and a 
non- viscous fluid 0. The HZ equation, on the other hand, arises in the study of massive nematic 
liquid crystals and in the study of shallow water waves In this paper we will show that the 
following system of two equations. 



^Jlu^^ - a- Xul^ 
and 

Uxxt = aux - 2uxUxx - uuxxx + ^ {ul)^^ , (4) 
are integrable for arbitrary values of the constant parameters a, A. Furthermore, they belong to 
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the same hierarchy corresponding to the negative and positive flows respectively, much hke the HD 
and the HZ equations. It is interesting to note from (0]) that when a 7^ 0, it can be scaled to unity 
through the scaling x x/a, u — > u/a. Therefore, looking at ^ alone, it would appear that 
there are only two meaningful values for a, namely, a = 0, 1. Similarly, under a scaling t t/a^^"^, 
X ^ x/a, u ^ u/a, the parameter a in the dHD equation (jS)) can be scaled to unity when it is not 
zero. However, there is no scaling which will scale a to unity (when it is nonzero) simultaneously 
in both the equations and, consequently, if the two systems belong to the same hierarchy, a has to 
be thought of as an arbitrary parameter. 

We note that, for a = A = 0, eqs. © and ((IJ) reduce respectively to the HD equation and 
the HZ equation Q. For lack of a better name, we will refer to © as the deformed Harry Dym 
equation (dHD) and Q) as the deformed Hunter-Zheng equation (dHZ). (The name "generalized 
Harry Dym" equation has already been used earlier in the context of a multi-component Harry 
Dym equation 9 .) For a ^ and A = equation @ was studied by Alber et al ^H] and its 
solutions contain solitons of the type known as umbilic solitons. Equation @ with a = 1 and 
A 7^ has appeared more recently in the literature TT| as describing short capillary-gravity waves 
and its solutions are known to become multi- valued in a flnite amount of time. While a matrix 
Lax pair for this system was provided in jllj and solubility of the model was argued based on a 
map to the sine-Gordon equation, interesting properties, such as the bi-Hamiltonian structure, the 
infinite number of conserved charges were not discussed at all. Our results clarify these aspects 
and provide direct support for the integrability of this new system. The deformed equation (jSJ for 
arbitrary a and A 7^ is truly a new integrable system (which to the best of our knowledge has not 
been studied in the literature) and leads to the above mentioned models in different limits. This 
can, therefore, be considered as a very rich system, much like the two boson hierarchy. 

Our paper is organized as follows. In section 2, we obtain the first Hamiltonian structure of the 
dHD and dHZ equations. In section 3, a second Hamiltonian structure, compatible with the first 
one, is obtained. We show that both the dHD and dHZ equations are bi-Hamiltonian and, therefore, 
integrable and belong to the same hierarchy corresponding to negative and positive flows. In section 
4, we present a Lax pair description of both the system of equations and construct the conserved 
charges of negative order from the Lax operator. We obtain a non-standard Lax description for 
the dHD equation for the special values of the deformation parameters A = 0, A = —3/ a. We 
argue that, for general values of the parameters, a non-standard Lax description will involve a 
pseudo-differential operator of inflnite order. In section 5 we present our conclusions. 
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2 dHZ and dHD as Hamiltonian Systems: 

To describe the dHD and the dHZ equations in a compact manner, let us introduce the following 
notation. Let us define 



F(c.x) = 2nxx -a- Xul^ , A = ^^^"^ . (5) 



Then, it follows from the definitions in ((SJ that 



K = {l-a\) = XF^ + (1 - Xu^^f = (A + A'^)F^ , 

_ 1 
F^ - 

-ps- = ■ (6) 

These and other relations following from these prove very useful in the analysis of the two systems. 
Given the dHZ equation @, we obtain from the definition in Q that 

A 



Ft 



F{u--ul 



(7) 



Similarly, in terms of the new variables in ((Sj), the dHD equation ^ can be written as 

ut = V2Ax , (8) 
and it follows that, under the evolution of dHD, 

Ft = V2AAxxx , 

At = -V2n^. (9) 

It is now straightforward to note from eqs. ((7| and Q that 

H_i = V2 j dxF (10) 

is conserved under both the dHD and dHZ flows. 

The dHZ equation can be obtained from a variational principle, 6 J dtdx C, with the Lagrangian 
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density 



1 0^2^4 / \ 

£ = -u^ut + + -uu^ - —u^ . (11) 



Ct 2 ^ 2 ^ A r, f 1 



(13) 



This is a first order Lagrangian density and, consequently, the Hamiltonian structure can be readily 
read out, or we can use, for example, Dirac's theory of constraints to obtain the Hamiltonian 
and the Hamiltonian operator associated with ()11|). The Lagrangian is degenerate and the primary 
constraint is obtained to be 

$ = vr - ^n, , (12) 
where vr = dC/dut is the canonical momentum. The total Hamiltonian can be written as 
Ht = j dx {irut - £ + /5^>) 
dx 

where /? is a Lagrange multiplier field. Using the canonical Poisson bracket relation 

{u{x),TT{y)} =6{x-y) , (14) 

with all others vanishing, it follows that the requirement of the primary constraint to be stationary 
under time evolution, 

mx),HT} = o, 

determines the Lagrange multiplier field /? in ()13() and the system has no further constraints. 
Using the canonical Poisson bracket relations ()14() . we can now calculate 

K{x, y) = {$(x), $(y)} = ]^dy5{x - y) - ]^dj{x - y) . (15) 

This shows that the constraint H12|) is second class and that the Dirac bracket between the basic 
variables has the form 



{u{x),u{y)}D = {u{x),u{y)} - j dz dz {u{x),^{z)}J{z, z'){^{z),u{y)} = J{x,y) 
where J is the inverse of the Poisson bracket of the constraint ((T3)) , 

J dzK{x,z)J{z,y) = 6{x - y) . 
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This last relation determines 



dxJ{x,y) = 5{x - y) , 



or 



J(x,y) = ViS^x - y) , 

where 

Vi = , (16) 

and can be thought of as the alternating step function in the coordinate space. We can now set 
the constraint l\12^ strongly to zero in (|13|) to obtain 

H2 = -Ht = j dx ^|n2 + ^un^ - ^u^^ . (17) 

Therefore, the dHZ equation can be written in the Hamiltonian form 

ou 

with Vi and H2 given by (fT^ and ((T7I) . respectively. 

Prom the results in we know that the HD and HZ equations belong to the same hierarchy 
of equations. Here, too, we will see that both the dHD and the dHZ equations belong to the same 
hierarchy. In particular, we note that 

5H_i 



Ut = Vi- 



6u 



with i^-i, given by (|1U() . yields the deformed Harry Dym equation (jSJ. As a result, the dHD 
equation also is Hamiltonian with the same Hamiltonian structure of the dHZ equation in ()16() 
and, consequently, has a Lagrangian description given by 

^ = ^UxUt - V2F . 

The reader can easily check that H2 is also conserved by both the dHD and dHZ equations, much 
like H^i. Note that for a = X = 0, these two charges reduce to the corresponding ones in the HD 
and HZ systems. 
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3 dHD and dHZ as bi-Hamiltonian Systems: 



It is well known that a system can be shown to be integrable if it is bi-Hamiltonian jl3l I14j . 
This corresponds to the system having a Hamiltonian description with two distinct Hamiltonian 
structures that are compatible. Therefore, we try to see if the dHD and the dHZ equations can be 
described as bi-Hamiltonian systems. For this, we have to find a second Hamiltonian description 
for the two systems. 

In order to determine a second Hamiltonian structure for the two systems, we recall (Tl that, 
when a = A = 0, the corresponding H^i is a Casimir of the second Hamiltonian structure for the 
HD and HZ systems, namely, 

^c=x=0) ^ Q^l ^ d-\xx (18) 

satisfies 

Since we know H^i for the deformed systems which is a generalization of H^^^^^^\ we look for a 
Hamiltonian structure for which it is a Casimir. With some work, it can be determined that H_i 
given in H10|l is a Casimir of 

V2 = Vt'^ = i (5-2f2 d-^ + d-^F^ 9-2) + A d-^uxx. d-^uxx.d-^ . (19) 

Note that this structure reduces to (jlHj) when a = A = 0. The skew symmetry of this Hamiltonian 
structure is manifest. The proof of the Jacobi identity for this structure as well as its compatibility 
with Vi in ()16() can be determined through the standard method of prolongation described in ref. 
jl4j . which we discuss briefly. 

Performing the change of variables 

W — Uxx 1 

the Hamiltonian structures (|16|1 and H19|l assume the forms 

V2 = \ [F^d + dF^) + \wx d~^wx . 
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We can construct the two bivectors associated with the two structures as 

e^i = ^ y dx {e^Vle} = ^jdxeA e^^^ , 

Using the prolongation relations, 

prux)je(u;^) = 2w pr vd^o^w) , 
prvx>^g{wx) = {prvx>^e{w))^ , 

prvx>2e{w) = —aOx + 'i.w 9x — XuP'Ox + WxO — XwWxO + \wx{d~^WxO) , 
prvD2e{w'^) = 2wprvD20{w) , 

prv-p^eiwx) = {prvv^eiw))^ , (20) 
it is straightforward to show that the prolongation of the bivector Qx>2 vanishes, 

pri;-p2e(6©2) = , 
implying that 1)2 satisfies Jacobi identity. Using H20|). it also follows that 

prVDj^g (6D2) + pr VT>29 (©Di) = , 

showing that Di and T>2 are compatible. Namely, not only are T>i,T>2 genuine Hamiltonian struc- 
tures, any arbitrary linear combination of them is as well. Any physical system that is Hamiltonian 
with respect to these two structures, therefore, defines a pencil system and is integrable. It is 
worth pointing out here that when a = A = 0, the second Hamiltonian structure (|18j) represents 
the centerless Virasoro algebra |15| (with dimension zero operators). The structure in ()19|) appears 
to be a highly nonlocal generalization of this algebra, but we are not familiar with any study of 
such an algebra in the literature. 

To show that the dHD and dHZ are bi-Hamiltonian systems, we note that the charges 

Hi = [ dxul, (21) 
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H-2 = Y^j dxFAl (22) 

are also conserved by both the dHZ and dHD equations. While the charge in (|2H) is the unmodified 
charge of the HD and HZ systems (this seems to be a simple coincidence) , the charge in (|22|) is a 
true generalization of the corresponding charge of the HD and HZ systems. With these charges, it 
is easy to see that the dHZ equation can be written in a truly bi-Hamiltonian form 

du ou 

Similarly, the dHD equation can also be written in the bi-Hamiltonian form 

ut = Vi— — = V2— — . 

du ou 

Thus, we see that both the dHD as well as dHZ systems are bi-Hamiltonian with the same two 
compatible Hamiltonian structures and are, therefore, integrable. 

4 The Lax Representation: 

When a system is bi-Hamiltonian, we can naturally define a hierarchy of commuting flows through 
the relation 

ni = ir„M=Pi^^=P2^, n = 0,1,2,.... (23) 

In the present case, both the Hamiltonian structures have Casimirs. We have already seen that 
H-i is a Casimir of P2 and it can be checked that the trivial Hamiltonian 



Ho = fd.u.^ (24) 



formally defines the Casimir of T>i (namely, if we write formally = 9^, it is annihilated by I?i). 
As a result, the system of flows can be extended to both positive and negative integer values for n. 
In this way, we see that much like in the HD and HZ systems [7j, the dHD and dHZ systems also 
belong to the same hierarchy corresponding to the negative and positive flows. 

Let us introduce the recursion operator following from the two Hamiltonian structures as 

R = V2V-^ . 
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Then, it follows from that 



and 



where 

B) = d'^Uxx + (-a + Uxx) <9~^ - A u^x d'^Uxx (26) 

is the adjoint of R. The conserved charges for the hierarchy can, of course, be determined in 
principle recursively from 1)25^ . However, in practice, integrating the recursion relation is highly 
nontrivial. Therefore, we look for a Lax representation for the system of dHD and dHZ equations 
which will allow us to construct the conserved charges directly. 

It is well known |16[ ll7j that for a bi-Hamiltonian system of evolution equations, ut = Kn[u], a 
natural Lax description 

is easily obtained where, we can identify 

M = , 

Here K'^ represents the Frechet derivative of Kn, defined by 

K'^[u]v = ^ Kn[u + ev] 
For the dHD and dHZ system of equations in ^ and respectively, we have 

i?^™ = K^2 = V2Kd^F-^d, 

B ^ = {-a + uxx) d-' + d(^u- ^ul^ . 
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The two systems have the same M = E) given in (|26() . It can now be checked that 

dM 



dt 
dM 



^ odHD 



^ „dHZ 



(27) 



dt 

do indeed generate the dHD and the dHZ equations and, thereby, provide a Lax pair for the system. 

One of the advantages of a Lax representation is that they directly give the conserved charges 
of the system. From the structure of (|27jl . it follows that TrM^*" are conserved, where "Tr" 
represents Adler's trace |18j. We note that 

2n+l 

TrM^- = , n > 1 
TrM5 = dxF , 



TrM"5 = — -I dxFAi , 
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TrM-t = - dx + -FAt 

n. \ F K 



(28) 



The first two nontrivial charges correspond respectively to H^i,H^2 given in eqs. (|lUj) and (|22|1 . 
constructed earlier by brute force. In fact, all H^n-i with positive n > can be constructed from 

2n—l ____ 

TrM ~ and by construction (namely, because of the nature of H27|) ). they are conserved under 
both the dHD and dHZ flows. Unfortunately, as is clear from H28p. this procedure does not yield the 
charges Hn with positive integer values. These are, in general, non-local and even in the HD and 
HZ case, construction of these charges relies primarily on the recursion relation (|25|) . It remains 
an interesting question to construct these charges in a more direct manner. 

The Harry Dym equation has a Gelfand-Dikii representation for the Lax pair, while the HZ 
equation does not. We will now discuss the existence of such a Lax representation for the dHD 
equation. A spectral problem associated with the dHD equations can be obtained from the recursion 
relation (|25p (see |19j and references therein) as follows. Introducing a spectral parameter and 
defining 



n=0 



n 

6u 
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we note that the recursion relation H25|) can be written compactly as (recall that Hq is a Casimir 
of Pi) 



(Pi - ;UP2)V'' = , 



or 



(1 - fiR^) -0^ = 



(29) 



which defines an eigenvalue problem for the eigenfunction ip'^ with eigenvalue A linear eigen- 
value problem can be derived from this if we can factorize the operator (1 — ^B)). 



Let us note that the operator B) in (|2(ij) can be rewritten in the form 



d-HFl,^+X)d-'+[Fl^,,+X)d 



-1 



3-2 



(30) 



where 



n=l 



and the coefficients X„ can be determined recursively to be 



^1 



0, 



When A = 0, it follows that X = and that F^^ ^^^^ is a simple function. In this case, the eigenvalue 
problem (|29|) can be f actor ized as 

(1 - i?t)^2 ^ 2 9-^-^503 [d' - f i^(Uo)) 'A = , 

where we have identified 
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This shows that if the hnear equation 







is satisfied, then will hold and this identifies the Lax operator for the system to be 

1 



L 



/?2 



In fact, it can be readily checked that when A = 0, the hierarchy of dHD equations can be obtained 
from the non-standard Lax equation 

dL 



dtn 



4\/2 



L,{L 



(2n-l)/2> 



>2 



The conserved quantities for this system can be obtained from Tr L(2"-i)/2, n = 0, 1, 2, . . .. 

On the other hand, when A 7^ 0, the coefficients Xn are nontrivial and X represents a pseudo- 
differential operator. The factorization, in such a case, is not so simple and, in principle would 
involve an infinite series of terms. For arbitrary values of k, the terms in the series can possibly 
be determined recursively. However, this is not very interesting. We simply note here that for the 
special value «; = 4, the infinite series of terms seems to have a simpler compact form and the Lax 
operator, in such a case, has the form 

F F 4 

and, in this case, the hierarchy of dHD equations can be obtained from the non-standard Lax 
representation 



dtn 



4\/2 



L 



(2n~l)/2 



K=4 



>2 



The conserved quantities, in this case, also follow from Tr L^^ 



(2n-l)/2 



and up to multiplicative factors. 



they have the forms given in (|28|) with k = 4. A simple Lax description for arbitrary A, however, 
remains an open question. 



5 Conclusion: 

In this paper, we have studied the general system of dHD and dHZ equations. We have shown that 
both these systems are bi-Hamiltonian and, therefore, integrable and belong to the same hierarchy 
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corresponding to negative and positive flows. The Lax pair for tlie two system of equations have 
been derived and conserved charges corresponding to negative integer values follow from the Lax 
operator. A simple construction of the charges for positive integer values remains an open question. 
For A = 0, we have constructed a non-standard Lax representation for the dHD equation, which 
involves a purely differential Lax operator. For arbitrary values of A, we have argued that a non- 
standard Lax representation will necessarily involve a Lax operator which is a pseudo-differential 
operator of infinite order. For the particular case of k = 4, however, this takes a simpler compact 
form. 
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